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■ The large- N c orientifold planar equivalence between TV = 1 SUSY Yang-Mills theory and 

ordinary 1-flavor QCD suggests that low-energy quark-gluon dynamics in QCD should be 
Cln' constrained by the supersymmetry of the parent theory. One SUSY relic expected from 



< 



orientifold equivalence is the approximate degeneracy of the scalar and pseudoscalar mesons 
in 1-flavor QCD. Here we study the role of the qq annihilation (hairpin) contributions to the 
meson correlators. These annihilation terms induce mass shifts of opposite sign in the scalar 
and pseudoscalar channels, making degeneracy plausible. Calculations of valence and hairpin 
correlators in quenched lattice QCD are consistent with approximate degeneracy, although 
the errors on the scalar hairpin are large. We also study the role of qq annihilation in the 1- 
and 2-flavor Nambu-Jona Lasinio model, where annihilation terms arise from the chiral field 
determinant representing the axial U(l) anomaly. Scalar-pseudoscalar degeneracy for the 1- 
flavor case reduces to a constraint on the relative size of the anomalous and non-anomalous 
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4-fermion couplings. 
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I. INTRODUCTION 

The emerging connections between QCD and string theory provide a promising new theoretical 
framework for the study of low energy hadron physics. A particularly interesting outgrowth of 
these developments is the large- N c orientifold planar equivalence between J\f = 1 SUSY Yang- 
Mills theory and ordinary nonsupersymmetric one-flavor QCD This correspondence can 
be described strictly as a field theory to field theory equivalence but its roots lie in string/gauge 
duality. To the extent that orientifold equivalence is reliable at N c = 3, it provides a powerful 
new source of insight into the chiral dynamics of quarks. What is particularly striking about this 
connection is that the quark of the daughter theory (one-flavor QCD) is the orientifold projection 
of a gluino in the parent theory (J\f = 1 SUSY YM). If this equivalence is even approximately valid 
at N c = 3, it would expose a deep and surprising role of supersymmetry in low energy quark-gluon 
dynamics. The prediction of orientifold equivalence for the size of the quark condensate in 1-flavor 



QCD has been compared with lattice results in Ref. [3]. In this paper we consider another of the 
predictions of orientifold equivalence - the degeneracy of the lowest lying scalar and pseudoscalar 
mesons of one-flavor QCD (hereafter referred to as a and if respectively) [4]- QCD itself provides 
no apparent reason to expect such degeneracy, while in the parent SUSY theory degeneracy is 
a simple result of the fact that the a and if belong to the same Wess-Zumino supermultiplet 
and that supersymmery is unbroken. Such a degeneracy, if valid in one-flavor QCD, would be an 
unmistakable relic of the supersymmetry of the parent theory. 

The properties of the scalar and pseudoscalar mesons of the parent theory are conveniently de- 
scribed in terms of the supersymmetric chiral Lagrangian constructed by Veneziano and Yankielow- 
icz (VY)[5]. In this formalism, the scalar and pseudoscalar mesons are associated with the real 
and imaginary parts of the complex scalar Wess-Zumino field representing gluino-antigluino bound 
states. From the viewpoint of QCD alone, the degeneracy of the scalar and pseudoscalar mesons 
seems somewhat mysterious. In that theory, we are accustomed to thinking in a large-iV c frame- 
work in which the pseudoscalar rf meson is a would-be Goldstone boson which acquires a mass via 
the axial anomaly. Lattice calculations have verified the role of the anomaly in generating the rf 
mass in detail by the study of the quark-line-disconnected "hairpin diagram" contribution to the 



pseudoscalar correlator 



(Here and elsewhere, we use the term "hairpin" to refer to a diagram 
involving quark-antiquark annihilation or creation where the quark and antiquark involved emerge 
from the same hadron.) In contrast to the would-be Goldstone boson interpretation of the r/, the 
scalar meson a, is, at least naively, a typical P-wave quark-antiquark meson. Phenomenology and 
lattice calculations suggest that the corresponding flavor nonsinglet mesons are quite heavy, ~ 1 
GeV. The flavor singlet a correlator includes hairpin diagrams along with the valence diagram, 
and, as we will discuss here, this results in a large negative mass shift of the flavor singlet relative 



to the nonsinglet. The scalar hairpin correlator was calculated in lattice QCD in Ref. 



as part 



of a study of the spin-parity structure of the OZI rule. In fact, the quenched calculations of Refs. 
[3, contain all the Monte Carlo results necessary for testing the orientifold prediction. Although 
errors on the scalar hairpin mass shift are sizeable, we will show that within these errors, the 
combination of the upward shift of the pseudoscalar mass and the downward shift of the scalar 
mass, when evaluated in the 1-flavor theory, brings the two mesons into approximate degeneracy. 

Througout this discussion, we adopt a large--/V c view of chiral dynamics, as typified by the 
Witten- Veneziano relation, which relates the hairpin-induced mass shift of the pseudoscalar meson 
to the topological susceptibility of quenched (pure-glue) QCD. In practical terms, we can take the 
derivation of the Witten- Veneziano relation as consisting of the assumption that the if correlator 
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in full QCD is well approximated by summing repeated quenched hairpin diagrams, with each 
hairpin vertex treated as a mass insertion. Although it would clearly be interesting to compare 
our results with the results of a full 1-fiavor QCD simulation, within the large N c framework the 
quenched hairpin diagrams actually provide a more direct measurement of the mass shift than the 
corresponding diagrams in full QCD, which measure a resummed combination of multiple mass 
insertions. For the pseudoscalar correlator, the interpretation of the hairpin vertex as a pure mass 
insertion has been confirmed in considerable detail by lattice calculations 

a a. 

Less detailed 

information is known for the scalar hairpin correlator, but we will assume that resummation of 
quenched hairpin vertices also provides a good description of the correlator in the scalar channel, 
so that the amputated, quenched hairpin correlator is the (mass) 2 insertion. It is interesting that, 
according to the Monte Carlo results [8] the scalar and pseudoscalar hairpin correlators are the 
only ones which are large, with the vector and axial-vector hairpins being more than an order 
of magnitude smaller. This reproduces the correct phenomenological pattern for the OZI rule, 
where, for example, the u-p mass splitting is tiny compared to the rj'-ir splitting. This entire 
OZI pattern can in some sense be seen as a relic of supersymmetry in that it reflects a similar 
structure in the Veneziano-Yankielowicz chiral Lagrangian for the parent SUSY theory. In the VY 
Lagrangian, the auxiliary complex field of the Wess-Zumino multiplet corresponds to F 2 + iFF in 
the gauge theory. In the low energy chiral Lagrangian, this field has no kinetic term, and, when 
integrated out, produces the mass terms for the scalar and pseudoscalar mesons, with the scalar 
mass arising from integration over F 2 and the pseudoscalar mass arising from integration over 
FF. Thus the gauge fluctuations represented by the auxiliary fields in the VY Lagrangian result 
in nonzero hairpin correlators only in the scalar and pseudoscalar channels, not in the vector and 
axial vector channels. This suggests that the complete pattern of observed OZI-rule violations 
in hadron phenomenology may be plausibly traced to the structure of the chiral Lagrangian in 
the supersymmetric parent theory. In particular, the importance of the scalar and pseudoscalar 
hairpin diagrams, and the relationship between them, reflects a fundamental SUSY-relic connection 
between the conformal anomaly and the chiral anomaly in QCD. 

Within the large- iV c framework, the quark dynamics leading to scalar-pseudoscalar degeneracy 
in 1-flavor QCD can be better understood by generalizing to the 2-flavor case. The separation 
between valence and hairpin terms is then obtained by considering appropriate combinations of the 
isotriplet and isosinglet correlators. In the 2-flavor theory we may identify the correlator given by 
the valence diagrams alone with that of the isotriplet meson, which in the pseudoscalar case would 
be a true Goldstone pion. Even in the 1-flavor theory, it is appropriate to think of the separation 
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between valence and hairpin graphs in the pseudo-scalar channel as separating the two types of 
chiral symmetry breaking which are taking place: the valence diagrams represent spontaneous 
xSB in the form of a Goldstone boson, while the hairpin diagrams represent the explicit breaking 
due to the chiral anomaly. In the VY Lagrangian, the pseudoscalar hairpin vertex is simply the rf 
mass term which arises from the integral over the auxiliary field FF. For the scalar channel, the 
separation between valence and hairpin contributions to the VY Lagrangian is not as obvious. In 
this case the valence diagrams also contribute to the scalar mass term, and in fact, the mass shift 
from the hairpin graphs turns out to be negative. 

In Section II, we analyze the Monte Carlo results obtained from scalar and pseudoscalar valence 
and hairpin correlators and discuss the significance of these results for the SUSY-relic prediction of 
scalar-pseudoscalar degeneracy in 1-flavor QCD. Although the evidence is far from conclusive, the 
qualitative pattern of valence and hairpin contributions to the a and rf masses which is required for 
degeneracy is clearly verified by the Monte Carlo results. We are thus encouraged to investigate 
the underlying quark interactions which lead to this pattern. In Section III we consider this 
question in a 2-flavor Nambu-Jona Lasinio (NJL) model. This model plausibly incorporates the 
chiral dynamics associated with short-range qq interactions. By studying the 2-flavor theory, we 
deconstruct the meson masses into their valence and hairpin components and show that the model 
predicts the same qualitative pattern of masses observed in the Monte Carlo results. In the 2-flavor 
NJL model, the axial U(l) anomaly is easily incorporated as an effective 4-fermion vertex which 
breaks the U{2) x U{2) symmetry of the NJL Lagrangian down to SU(2) x SU(2) Qj. The mass 
shifts resulting from the axial anomaly can thus be easily determined by a simple application of 
the standard Bethe-Salpeter analysis. We then discuss the effective Lagrangian and mass spectrum 
of the 2-flavor theory. The effective meson Lagrangian contains two ^] chiral field multiplets 
(under SU{2) x SU(2)), describing both 1 = and 1=1 mesons. The mass terms in this theory 
provide the information needed to show that the non-anomalous U(2) x U (2) invariant 4-fermion 
interaction C s represents pure t-channel glue exchange between valence quarks, while the term C a 
representing the anomaly is a pure hairpin or qq annihilation term. Within the large- N c framework, 
once we have identified the valence-hairpin separation we can generalize the 2-flavor results for the 
correlators to any number of flavors by including a factor of Nj/2 for each annihilation vertex. 
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FIG. 1: Scalar and pseudoscalar mass 2 for one-flavor QCD (MeV 2 ) 
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Total Mass 2 
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+ [407(11)] 2 


+ [515+ 12- 13] 2 


0++ 


+ [I4I6(I4)] 2 


-[1350(90)] 2 


+ [427 + 249 - 756] 2 



TABLE I: One-flavor QCD results for scalar and pseudoscalar meson masses. The columns list values in 
Mev 2 for m Vl mh v , and m, respectively. Errors quoted are statistical only. 

II. MONTE CARLO RESULTS 



Ideally one would like to test orientifold equivalence in the chiral limit where it does not suffer 

n 

from finite quark mass corrections. Since the accuracy of the data for the scalar hairpin [8] does 
not permit a reliable chiral extrapolation, we test the prediction of scalar-pseudoscalar degeneracy 
using one of the lightest quark masses employed in the study of Ref. [8j (m n = 315(6) MeV), 
corresponding to clover improved Wilson fermions with C sw = 1.57, hopping parameter k = .1427, 
using 300 quenched gauge configurations generated with Wilson action at (3 = 5.7 on a 12 3 x 24 
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lattice. 

The valence or nonsinglet mass was extracted from combined exponential fits to smeared and 
local scalar and pseudoscalar valence correlators, 

f 2 

d 3 x{u(x)r i d(x)d(0)r i u(0)) ~ Jj^e-m^xo /-q 

2m v 

where u and d are two different flavor Dirac spinors, and Tj = 1 or 75. For the scalar valence 
correlator a quenched artifact associated with the v'-tt intermediate state must be removed from the 
correlator before performing the exponential fit [7]. The result for m ao is in reasonable agreement 



with recent full QCD calculations 



id]. 



The mass shift in the singlet channel is extracted from the hairpin correlator, 

^hp(x) = (u(x)Tiu(x)d(Q)Tid(Q)) c (2) 

Here the subscript c denotes the connected correlator, i.e. for Tj = 1 a constant (VEV) 2 must 
be subtracted off. The Monte Carlo results for the quenched hairpin correlators are well fit by a 
double pole formula 



/ 



2™2 



f . 777, 

d 3 xA hp (x) ~ + m„x )e-~ (3) 



Upon resummation of multiple hairpin insertions, the mass of the flavor singlet meson is shifted 
from the nonsinglet value m v to m, where 

m 2 = m 2 v + m 2 hp (4) 

The results of these studies for k = .1427 are summarized in Table U and Fig. 1. The most 
striking finding is the large negative contribution of the hairpin to the mass of the scalar meson, as 
shown in Fig. 1. Although the statistical errors in the measurement of the scalar hairpin diagram 
are large, the size and sign of the scalar mass shift make scalar-pseudoscalar mass degeneracy 
plausible. 



III. CHIRAL LAGRANGIANS IN ORDINARY AND SUSY QCD 

The pattern of chiral symmetry breaking in J\f = 1 SUSY Yang-Mills theory with an SU(N) 
gauge group is well understood. The usual description starts with the chiral anomaly breaking 
the U(l) chiral symmetry of the gluinos down to Zin- The discrete Z^n chiral symmetry is then 
spontaneously broken down to Z2 by the appearance of N possible discrete gluino condensates 
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in the vacuum. In this description, one obtains domain walls between chiral condensates, but no 
Goldstone boson. Nevertheless, one may assume that, for large N, the fluctuation of these domain 
walls can be approximately described by a continuously varying order parameter. Adopting this 
viewpoint, the Veneziano-Yankielowicz Lagrangian generalizes Witten's interpretation of the rf 
meson in QCD as a would-be Goldstone boson [11:] to the case of J\f = 1 SUSY YM. For the 
SUSY case, the pseudoscalar meson must then be accompanied in the chiral Lagrangian by the 
other members of its supermultiplet, which are degenerate with it since supersymmetry remains 
unbroken, even in the presence of a chiral condensate. (The broken chiral U(l) symmetry is R- 
symmetry, which is not a supersymmetry transformation, but rather an outer automorphism of 
the SUSY algebra.) In addition to the complex scalar field cj>, representing the gluino condensate 
AA, the basic supermultiplet of the chiral Lagrangian includes a fermion field representing gluino- 
antigluon bound states and a complex scalar auxiliary field which is identified with F 2 + iFF of 
the gauge field, where 

F 2 = F- u F a ^, FF^-e^F^ (5) 

The fact that the gauge densities F 2 and FF appear in the Lagrangian as auxiliary fields with 
no kinetic term amounts to a "heavy glueball" approximation in which effects associated with real 
propagation of glueballs are neglected. The auxiliary fields in the chiral Lagrangian reflect the 
relevant vacuum structure of the gauge field. 

The construction of the chiral Lagrangian for J\f = 1 SUSY YM is discussed in detail in the 
original reference [3]. Here we summarize the terms of interest in our discussion. The WZ chiral 
field is given in terms of the gluino condensate order parameter by 

4> = cXr\ l , <t>* = cX L \ R , (6) 

Since we are only discussing the bosonic sector, we will not need the fermionic superpartner field 
which is a gluino-gluon composite. We consider only the terms in the Lagrangian that refer to the 
field (f) and the auxiliary field 

M = — ^c(F 2 + iFF) (7) 
In ([6]) and (JT]), the constant c is given in terms of the renormalization group beta function 

c = (3(g)/2g = -^9 2 N c /tt 2 + 0(g A ) (8) 
The auxiliary fields are coupled to <j> so that, when they are integrated out, it sets them equal to 

cF 2 = -^a(<^) 2/3 log (0*0///) + • . . (9) 
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and 

cFF = Ua(<f)*<j)) 2/3 log (070) + . . . (10) 
where . . . represents terms that we are not considering here. The resulting potential for is 

V S M = ^a(0*0) 2 / 3 log(0//i 3 ) log(07^ 3 ) (11) 
This potential generalizes to the supersymmetric case the large N chiral Lagrangian formulation 



111 ] , in which the chiral anomaly is represented by an rj mass term constructed from the log of the 



chiral field determinant. An exponential parametrization 

= M 3 e (<T+ir/)// (12) 

(with / = pseudoscalar decay constant) shows that the potential (jlip produces a supersymmetric 
mass term for the scalar and pseudoscalar mesons, 

Vss = + r/ 2 ) + ... = \m*„{p* + r/ 2 ) + . . . (13) 

IV. VALENCE AND HAIRPIN MASSES IN THE NAMBU-JONA LASINIO MODEL 

As discussed in Section II, the lattice determination of the flavor-singlet meson masses requires 
two separate calculations for the valence and hairpin terms in the correlator. In the pseudoscalar 
channel, the valence diagram is a massless Goldstone boson propagator, and the hairpin vertex 
is a positive mass insertion. By contrast, the scalar valence correlator gives a large nonsinglet 
mass, and the scalar hairpin is a negative mass insertion. In this Section, we will investigate the 
various contributions to the scalar and pseudoscalar masses in the Nambu-Jona Lasinio (NJL) 
model. This model has proven to be a useful approximation to QCD which retains much of the 



chiral structure of the full theory. T 



application to QCD phenomenology [ 



rere are several modern reviews of the NJL model and its 
\\\ . Here we only briefly summarize properties of the 



0, Q,C 



model relevant to determining scalar and pseudoscalar masses. In order to sort out the valence 
and hairpin contributions in the chiral Lagrangian of 1-flavor QCD, it is useful to generalize to the 
2-flavor theory. We define an / = I doublet of quark spinors, 



(14) 
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We then write the basic 4-quark interaction of the NJL model in terms of the following quark 
bilinears: 

£ = qq (15) 

= qi^nq (16) 

N = qi 75 q (17) 

A* = qnq (18) 

The four-vectors 

5=(S,n) (19) 

and 

P = (N, A) (20) 

each transform as a [i, |] representation of chiral St/ (2) x SU{2). An invariant 4-fermion interaction 
is given by an arbitrary linear combination of the two invariants, 

A = s 2 + n 2 (21) 
C 2 = N 2 + A 2 (22) 

The U(2) x U{2) invariant NJL interaction is given by 

C s =g s {C 1 +C 2 ) (23) 

Like the QCD lagrangian, C s is invariant under an axial U(l) transformation q —> e ia,ys ' 2 q, or 

(S±iP) -» e ±ia (S±iP) (24) 

The axial t/(l) anomaly can be represented by an effective 4-fermion interaction which is invariant 
under SU(2) x SU(2) but violates the 17(1) symmetry (f2"4"j) [9]. Defining the 2x2 matrix chiral 
field 

$ab = q a {l ~ lb)qb (25) 

we take 

£a = 2g a (Det<£> + h.c.) = g a {C\ - £ 2 ) (26) 
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FIG. 2: Quark sclf-cnergy 





+ 



FIG. 3: Resummed bubble diagrams for the meson correlator 



+ ... 



Note that the chiral field can also be written as 



(27) 



Results of the NJL model are given in terms of a constituent quark mass m q , which is the 
solution to a self-consistent BCS-type gap equation for the quark self-energy depicted in Fig. 2. 
The mass of the meson in a given spin-parity channel is obtained from the solution to a Bethe- 
Salpeter equation, which amounts to a resummed bubble approximation for the meson correlator 
as shown in Fig 3. A bound state pole corresponds to a zero of the inverse correlator. Let us first 
consider the meson masses that are obtained if we neglect the axial anomaly term and include only 
the C s interaction. Then the gap equation is 



1 = g s Iv 



(28) 



where 



L, 



JAmi 



1 




(29) 



2vr 2 j 4m 2 

The meson correlators are given by a resummation of the one-loop integral in Fig. 3, which reduces 
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to a dispersion integral of the form 



i 



F (M^r^i-^f* (3„) 

2tH J 4m 2 s-q 2 y s J 
Let us denote the scalar meson correlators by D a for the 1 = channel and D ao for the 1=1 
channel, and the pseudoscalar meson correlators by for the 1 = channel and D w for the 1=1 
channel. In the theory without the Ua(1) anomaly, the inverse correlators are (dropping an overall 
constant) 

D- / 1 (q 2 )=D- 1 ( q 2 ) = l- gs F(q 2 ,0) (31) 

and 

D-\q 2 ) = D^(q 2 ) = 1 - g s F(q 2 ,4m 2 q ) (32) 

Noting that F(fi 2 , ^i 2 ) = I vac for any // 2 , the zeroes of the inverse correlators give the meson masses, 

m 2 , = m 2 = (33) 

m l = m l = 4m l ( 34 ) 

The inclusion of the axial anomaly term C a simply changes the relative strength of the 4-fermion 
couplings in the various channels, but the Bethe-Salpeter bubble resummation analysis of C s + C a 
is similar to the previous discussion. The mass gap equation becomes 

1 = {9s + 9a)Ivac (35) 

Here g a = (1 + ^)g a includes a Fierz-transformed cross contraction term. The inverse correlators 
for the four channels are now 

Zr,V) = l-(g s -~g a )F(q 2 ,0) (36) 

D-\q 2 ) = l-(g s+ ~g a ) F (q 2 ,0) (37) 

D-\q 2 ) = l-(g s + ~g a )F(q 2 ,4m 2 ) (38) 

^ooV) = 1 " (9s ~ ~g a )F{q 2 Am 2 ) (39) 

Note that the coupling constant for the isovector ir channel is the same as it is for the gap equation 
(j35|) . so again using F(0, 0) = I vac we see that the pion remains massless, as it should. The a 
channel also has the same coupling as the gap equation, so the mass of the a remains as before 
when expressed in terms of the quark mass, 

ml = 4m 2 (40) 
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For the and ao masses, the difference in coupling from the gap equation results in an additional 
mass shift induced by the Ua(X) anomaly. It is convenient to write the associated Bethe-Salpeter 
equations in terms of the function 



i 

,2\ 2 



/(9 2 ) = A/ -^(l-^l ds (41) 
4vr 2 J 4m 2 s-q 2 y s J 

The resulting equations for the masses of the rj' and ao mesons then reduce to 

K'/K) = (42) 

and 

(<-^f(ml) = ±-J^ (43) 



These equations can be easily solved numerically. For the present discussion, it is sufficient to 
consider the approximation in which the cutoff mass A 2 is large compared to q 2 and m 2 , in which 
case 

/«) « /«) « /(0) « In ^ + const. (44) 

In this approximation, the meson masses are 

ml = 0, m 2 , = ml (45) 
m 2 a = 4m 2 , m 2 = 4m 2 + m\ (46) 

where we have denoted by m\ the mass shift associated with the chiral £7(1) anomaly, 

2 1 ( 2 <7a 



'"-^tmW^iJ (47) 

Although it appears that the effect of the U (1) anomaly in the scalar sector is to leave m a unchanged 
and increase m ao , the interpretation of (|46p must take into account the fact that the anomalous 
g a term also enters the mass gap equation and redefines the quark mass. As shown below by a 
bosonization argument, the actual effect of the U{\) anomaly term in the scalar channel is to leave 
the ao mass unchanged and lower the a mass. This identifies C a as an annihilation term which 
only affects the isosinglet channels. 

The pattern of scalar and pseudoscalar masses which emerges from the NJL model can be 
understood more generally by considering the effective meson Lagrangian. To proceed, we need 
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to replace the quark bilinears (fT5j) by meson fields a, jf, rf, and Sq. This can be accomplished by 
writing an exponential parametrization of the chiral field <1>, Eq. ([2] 



$ = /i 3 exp[(cr + in' + do • r + in ■ f)/f] (48) 

If we ignore the Z7a(1) anomaly, the ?7(2) x f7(2) invariant model can be described by an effective 
potential of the linear sigma model form, 

V = — %Tr&$> + -^Tr(&$) 2 (49) 
/i 2 2fjfi 

Expanding out to quadratic terms in the meson fields, one finds equal mass terms for the a and ao 
fields, 

V « const. + ^ml(cr 2 + Oq) (50) 

where m 2 = 4a/x 4 // 2 . Thus we get the same mass spectrum predicted by the Bethe-Salpeter 
analysis of the anomaly-free NJL model, Eqs. (|33|)- (|34|) . The degeneracy between isosinglet and 
isotriplet mesons in each channel (m a = m aQ and rry = m n = 0) indicates that, in the absence of 
the axial anomaly, the U(2) x U(2) invariant NJL interaction (|23|) includes no hairpin annihilation 
contribution. The meson propagators obtained in this case are the analog of pure valence prop- 
agators in QCD. The absence of hairpin contributions in C s is a specific feature of the 4-fermion 
interaction, and not a consequence of U{2) x U{2) invariance. To see this, note that in the chiral 
Lagrangian (|49p . we could have included an invariant term of the form 

(Tr&§) 2 (51) 

Although it is U(2) x U(2) invariant, such a term breaks the degeneracy between a and ao masses, 
and thus represents a hairpin annihilation term in the QCD scalar meson propagator. 

To include the axial U(l) anomaly into the effective meson Lagrangian for the 2-flavor theory, 
let us first consider the bosonized form of the determinant interaction (I26p . Using (|48j) we find 

g a [Det<& + h.c] = g a a 6 [ e 2 (-+^')// + e 2(<-Wj ^ 2g y ^ + 2a/f + 2(fJ 2 _ ^yf + ^ (52) 
Thus, the chiral anomaly introduces a negative mass shift for the a along with the expected positive 



rj' mass. 



V^V + ^m 2 ( V ' 2 -a 2 ) (53) 
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where uiq = 8g a/ u 6 // 2 . This reproduces the same pattern of masses (|45p -([46p found in both the 
lattice Monte Carlo calculations and in the NJL model with the anomaly term included, 



M 



2 



0, 



m 2 , 



m 



Cl() 



M 2 



(54) 
(55) 



The bosonic formulation of the effective Lagrangian makes it more obvious that the effect of the 
chiral anomaly in the scalar sector is a negative shift of the a mass, not a positive shift of the 
ao mass. In the NJL bound-state analysis, this requires an interplay between the ao bound state 
equation and the gap equation defining the quark mass. 

There is a well-known shortcoming of the form of the chiral anomaly represented by the deter- 
minant (|26p \\\ . In its bosonized form, the determinant contains higher order terms in 
the a and rj' fields corresponding to multiply OZI-suppressed meson emission vertices which should 
be down by additional powers of 1/N C . The form of the anomaly that satisfies proper OZI counting 
for large N c is a pure mass term constructed from the log of the determinant, 



(logDet($/^ 3 )) 2 + h.c 



1 



-ml(r]'' 



a 



(56) 



Comparing this with (I53p we see that the large- N c form of the chiral anomaly, with the appropriate 
coefficient, produces the same mass shifts as the 4-fermion Det<& form used in the NJL analysis. 



V. DISCUSSION 



From the results of the 2-flavor NJL model, it is straightforward to separate the valence from the 
hairpin contributions to the scalar and pseudoscalar meson correlators. The U(2) x U(2) invariant 
interaction C s , Eq. (|23|) . represents a pure valence interaction (i.e. the qq interaction arising from 
t-channel glue exchange) with no qq annihilation. The hairpin annihilation vertex can be identified 
with the anomaly term £ OJ Eq. ([25]). which breaks U{2) x U{2) down to SU{2) x SU{2). This 
separation is inferred from the fact that, if only C s is included in the Bethe-Salpeter analysis, there 
is no splitting between the 1 = and 1 = 1 mesons in either channel, identifying this term as 
a pure valence interaction. When the C a interaction is included, it effects only the flavor-singlet 
channels, identifying this term as a pure qq annihilation. (Note that in the NJL approximation, 
even if the actual qq interaction is a t-channel exchange valence graph, a Fierz transformation of 
the 4-fermion vertex allows us to write the spin contractions in the s-channel form, so that the 
spinor products factorize into traces around each loop in the Bethe-Salpeter bubble sum.) 
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As shown in Section IV, in the meson Lagrangian, the effect of the anomaly interaction C a in 
the scalar channel is to introduce a negative mass shift for the a. This agrees at least qualitatively 
with the Monte Carlo results for the hairpin mass shifts. Quantitatively, the chiral Lagrangian 
analysis predicts equal and opposite (mass) 2 shifts, Eq. ([56]) . while the Monte Carlo results in 
Table I and Fig. 1 suggest a much larger shift in the scalar channel. However, with the large errors 
in the scalar Monte Carlo data, the significance of this discrepency is not clear. A more accurate 
study of the scalar hairpin correlator would clearly be of great interest. 

One of the broad goals of the present study is to expose the details of the quark dynamics 
which underly the SUSY relic prediction of scalar-pseudoscalar degeneracy. [In this regard, it is 



interesting that recent work on the holographic brane construction of QCD from string theory 
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has derived a generalized NJL model for quark dynamics starting with a construction of chiral 
quarks based on D8 branes It is instructive to compare the effective Lagrangian treatment 
of the NJL model in Section IV with the supersymmetric Veneziano-Yankielowicz Lagrangian 
which is presumed by orientifold equivalence to approximately describe 1-flavor QCD. In the VY 
Lagrangian, a-rj degeneracy is a consequence of the fact that the conformal anomaly oc F 2 and 
the chiral anomaly oc FF are in the same supermultiplet and hence have equal couplings to the 
log (j) terms (c.f. ([9]) and (jlOp ). Thus in QCD, the statement of scalar-pseudoscalar degeneracy may 
be interpreted as a relation between the relative strength of the valence and hairpin interactions. 
Applying the analysis of Section IV to the 1-flavor theory gives an equation for my. The scalar 
mass is given by the usual to 2 = Am 2 . So if we impose the SUSY constraint 

ml = m 2 , (57) 

we find that a-rj' degeneracy implies that the g s and g a couplings are related by 

f x (2s - l)dx „ f x dx 

9s / / 7 -in = 9a / / i T\ ( 58 ) 
Jl y x(x — 1) Ji y x{x — 1) 

Where g a = g a /2 = 3g a /4: and A = A 2 /4m 2 . This constraint determines the strength g a of the 
anomaly term needed to give a-rj' degeneracy. We note that, for any cutoff A 2 > m 2 ,, Eq. (|58p 
gives g a < g s . For a cutoff close to threshold, A ~ mo, it gives g a ~ g s , while for large cutoff 
A >> to , we find g a « g s . 

The valence-hairpin separation realized in the 2-flavor theory also applies to the 1-flavor the- 
ory, and we can identify the U(l) invariant and U(l) violating 4-fermion interactions as valence 
and hairpin terms, respectively. Moreover, we can relate this separation to a separation between 
integrals over the F 2 and FF auxiliary fields in the VY Lagrangian. From ([9]), (jlOp . and (|12p . it is 
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seen that the sigma mass comes entirely from the F 2 term and the rf mass entirely from the FF 
term. Thus the FF integration contributes only to the pseudoscalar hairpin interaction, while the 
F 2 integration contributes to both the valence and hairpin interactions in the scalar channel. We 
can identify the correct valence-hairpin separation in the 1-flavor VY effective potential (llip . 

V ss = Vval + Vhp (59) 

The hairpin term can be identified with the large- N c form of the chiral anomaly (|56p . which for 
the 1-flavor case is given in terms of the £7(1) chiral field (]12p by 

VhM = ^a(<f 0) 2/3 [ilog(0//z 3 )] 2 + h.c. (60) 

Subtracting this from the supersymmetric potential, we see that the valence term in the VY 
Lagrangian is given by 

V va i(4>) = ^«(0» 2 / 3 [log(«f <£/M 6 )] 2 (61) 

This gives the correct separation of mass terms, with V va i = m^a 2 and Vhp = ^ m o(v' 2 ~ (j2 )- 

This work was supported in part by the Department of Energy under grant DE-FG02- 
97ER41027. 
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